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Abstract 

We apply the time-dependent local density approximation (TDLDA) to cal- 
culate dipole excitations in small carbon clusters. A strong low-frequency 
mode is found which agrees well with observation for clusters C„ with n in 
the range 7 — 15. The size dependence of the mode may be understood sim- 
ply as the classical resonance of electrons in a conducting needle. For a ring 
geometry, the lowest collective mode occurs at about twice the frequency of 
the collective mode in the linear chain, and this may also be understood in 
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I. INTRODUCTION 



In this work we examine strong optical transitions in light carbon clusters, with the goal 
of providing a diagnostic tool for the structure of the clusters. Carbon has a rather complex 
structural evolution as a function of the number of atoms n. Small clusters are linear, 
and at higher n rings are favored. For intermediate n values, rings may be more favored at 
even n and linear chains at odd n. Direct information on the shapes of the clusters may be 
obtained by ion chromatography 0], but one mostly relies on spectroscopic data to interpret 
the structure [^]. 

We will demonstrate that the strong optical absoption resonances behave in a very sys- 
tematic way as function of cluster size, and that the frequency of the transitions gives direct 
information about the shape of the cluster. In particular, the strong tt — tt* transition in 
chains is at about half the frequency of the mode in rings. We also find that the frequencies 
can be understood very simply using classical concepts. The most important determinant 
in systems of similar composition is the polarizability, which is much larger for chains along 
their axes than for rings in the plane of the ring. 

In the section 2 below we discuss our quantum calculations, which are done using the 
time-dependent local-density approximation (TDLDA). This is nothing more than the Kohn- 
Sham equations with the "energy" of the orbital replaced by the time derivative id/dt. The 
TDLDA and its small amplitude limit has been shown to be a reliable tool to calculate 
strong resonances in the optical response of metallic clusters and molecules such as 
polyenes [0]. Our method to solve the equations is a very straightforward one, but it deserves 
some discussion because there are many implementations of small-amplitude TDLDA in the 
literature, and our method is not conventional in condensed-matter or chemical physics. 
We use a Hamiltonian of the usual Kohn-Sham form, taking a local-density approximation 
for the non-Coulombic electron-electron interaction, and a pseudopotential to describe the 
carbon ions. 

However, the calculations with ah initio Hamiltonians are very time consuming, and it is 
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important to analyze the results to gain understanding of the functional dependence of the 
excitation energy on the cluster properties. This will be done in Sect. 3; our result briefly 
is that the resonance frequency Un in chains scales with n as oc yhin/n. This applies to 
both chains and rings, with different coefficients of proportionality. 

Part of the interest in the specta of carbon chains comes from astrophysics, the question 
of the composition of interstellar matter. Absorption bands are seen which may be due to 
carbon clusters, but specific identification of species remains controversial . The TDLDA 
is probably only reliable to 10% or so on the frequencies of the strong transitions, which is 
not enough to use the theory to assign unknown transitions. However, it might be that the 
theory could be used to predict the shifts due to small changes in structure, for example by 
adding another atom at the end of the chain. We plan in the future to investigate various 
molecules that are similar to the carbon chains and see how large the perturbations are. 

II. TDLDA 
A. General 

Time-dependent mean field theory is a powerful tool to calculate the collective excita- 
tions of quantum many-particle systems, and it has been widely applied in cluster physics. 
There are a number of implementations of the theory, all starting from the time-dependent 
Schroedinger equation. With the local density approximation for the electron-electron in- 
teractions, the equation has the form of a time- dependent Kohn-Sham equation, 

zh^^Ur,t) = ||^V2+ E VUr-Ra) + e'J rfr"^^^ + V;,(p(r^ t))U,(f, t). (1) 

Here ipi represents a single-particle wave function, n is the electron density, V^on is the ionic 
core potential, and V^c is the potential associated with exchange and correlations. For our 
purposes, ipi is not perturbed greatly from the ground state wave function, and a small 
amplitude approximation justifiable. This leads to the RPA and the linear response method 
of solution. We believe that for large numbers of particles, the most efficient technique to 
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treat LDA Hamiltonians is the straightforward integration of the time-dependent equations 
of motion. The argument is given in the Appendix. The optical response of the hght chains 
C3.5 and C7 has been calculated with the configuration-interaction method of quantum 
chemistry [|1^,|11[], but this brute- force technique is impractical in large clusters. 

Our numerical method is taken over from the nuclear physics []T^, where it is called 
the TDHF theory. The efficiency of the method requires that the Hamiltonian matrix be 
sparse, and this is fulfilled by using a coordinate mesh to represent the wave function. It 
is also necessary for numerical efficiency to use Hamiltonians whose energy scales are not 
too large; this means that one treats only valence electrons and uses pseudopotentials for 
Vion to take into account the effects of the core electrons. Thus we treat dynamically only 
the four valence electrons of carbon. Our psuedopotentials are calculated according to the 
commonly accepted prescriptions ||13|Jl^ . The exchange-correlation potential Vxc is taken 
from ref. ||15|| . 

Some details of our integration algorithm are given in the Appendix. The important 
numerical parameters are the spatial mesh size, Ax, the number of mesh points M, the time 
step At, and the total length of time integration T. The values for these parameters to C„ 
clusters are discussed below. 



B. Initial conditions 

We wish to start with the Kohn-Sham ground state wave function, and in principle the 
geometry of the ions as well as the electron wave functions are determined by minimizing 
the Hamiltonian function. However, a full minimization of the structures is quite time 
consuming, and we believe that small variations in bond lengths will not have a significant 
effect on the collective excitations. We therefore take geometries from outside rather than 
calculating ab initio. For simplicity, we fixed the nearest neighbor distance of the carbon 
atoms at 1.28 A. This is the average LDA equilibrium distance for large circular rings or 
long chains. 
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To calculate the response, our initial wave functions are perturbed from the static solu- 
tions by a velocity field 



Uf,0)^e^'^ct>i{f^. (2) 

The momentum k is set to be small to ensure that the response is linear (we use typically 
k = 0.1 A~^). The real-time evolution of the dipole moment is obtained as 

z{t)=J2<iP,{t)\z\i^^it) >, (3) 
i 

and its Fourier transform in time gives the dipole strength function S{uj) 

S(u)^ — — / dtz(t)sm(ujt). (4) 
Trk Jo 

The strength function defined in this way is related to the oscillator strength / by 

I duS{u) = /. (5) 



C. Results 

We first discuss the numerical parameters needed for carbon clusters. We found that a 
mesh size of Aa; = 0.3 A is needed to adequately represent the sp-orbitals of carbon. This 
is more than a factor of two finer than the mesh size parameter needed for the s-orbitals 
of alkali metals. The wave function is represented on this grid within the interior of a 
cylindrical volume. The size of the cylinder is chosen to include all mesh points within 4 
A of each carbon atom. With this procedure, the wave function for the Cy cluster requires 
30,000 mesh points. 

There are two numerical time parameters. The first is the time step in the integration. 
In the carbon calculations we used At = 0.001 eV~^, which is an order of magnitude smaller 
than the value needed for alkah metal clusters. This corresponds to the order-of-magnitude 
difference between the maximum kinetic energies in the meshes used for systems, 0.3 A in 
carbon versus 0.8 A in alkali metals. The other numerical time parameter is the total 
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integration time T. The inverse of this sets the scale for the energy resolution in the response; 
we use T = 10 eV~^ to make visible structure on the energy scale of 0.1 eV. 

The response of a typical case, C7, is shown in Fig. 1. Plotted is the integrated oscillator 
strength as a function of excitation energy. The upper curve shows the response parallel 
to the chain, and the lower curve shows the perpendicular response. The total sum rule 
(/ = 4n = 28, four valence electrons per carbon atom), is satisfied to within 10% by the 
calculation^ The longitudinal response has a strong excitation at 5.3 eV. About 1/3 of the 
28 units of oscillator strength is in this excitation. The perpendicular response has no 
corresponding strong excitation below 10 eV. The small wiggles in the calculated response 
are artifacts of the truncation in the Fourier transformation. In our case here, the effect 
is obviously small. To understand better the nature of the 5.3 eV collective excitation, 
it is helpful to know the characteristics of the orbitals near the Fermi surface and the 
corresponding single-particle transitions. The single-particle level scheme, giving by the 
eigenenergies of the static Hamiltonian, are shown in Fig. 2. The lowest strong single- 
particle transition is between the highest occupied and lowest empty tt orbitals. This gap is 
2.1 eV, roughly half the frequency of the TDLDA mode. Note that there are two a orbitals 
within the vr — vr* gap, but these electrons are localized on the outer ends of the carbon chain 
and do not couple strongly to the other states. 

The 2.1 eV tt — vr* transition has a very large oscillator strength (/ = 10.2), which equals 
the strength of the TDLDA mode to within 4%. The magnitude of / may be understood 
qualitatively as equal to the number of vr electrons. To a good approximation, the longitu- 
dinal dipole operator on vr electrons in a chain keeps them in the vr manifold, so the sum 
rule is approximately conserved within the vr manifold. The number of them in an odd-n 



^In principle, the TDLDA conserves the sum rule exactly if the Hamiltonian is local. This is 
not the case when the core electrons are omitted, because the resulting pseudopotentials have a 
nonlocal character [14|. 
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chain is = 2n — 2, filling the molecular orbitals in the usual way. Thus, if the lowest tran- 
sition would exhaust the sum rule, the oscillator strength would be / = 12 for the twelve ir 
electrons in Cj. This is to be compared with ~ 10 from the LDA wave functions. 

We have calculated the excitation energies and oscillator strengths for the collective 
transition for carbon chains in the rangeQ = 3 — 20, and the results are tabulated in Table 
I. When more than one transition is seen, the average excitation energy and the summed 
oscillator strength are given. The features described for Cj apply systematically these chains. 
The oscillator strength is roughly given by the number of vr electrons; the frequency of the 
TDLDA mode is larger than the unperturbed excitation energy by a factor that ranges from 
2 in the light chains to 3 in the heavier ones. The frequencies from Table I are graphed in 
Fig. 3 and one can see that they vary smoothly with chain length. The filled circles are 
experimental points for odd n from ref. |jl6[ and [|l^. We see rather good agreement for 
the larger chains. The clusters Ca^sj have been calculated by the CI method of quantum 
chemistry [|11|, and the predicted strong transition is in the same region, at 8.10, 6.35, and 
5.54 eV, respectively, which are close to our results of 8.1, 6.4, and 5.3 eV. The oscillator 



strength for C3 is calculated in ref. W^ as about 1.1. The oscillator strength given in Table 
I is 3.1 for C3, but this is the strength along the chain direction. Dividing the strength 
by three to average over spatial orientations gives f=1.0, in agreement with the quantum 
chemistry calculation. While there is experimental data for the excitation energies of odd-n 
chains in the range n = 7 — 15, there is no corresponding data on the even chains. We are 
confident that the systematics is smooth going over odd and even n, so it should be possible 
to observe the mode in even-n. 



Ref. [17] also reported weaker transitions in odd-n chains at lower energy. The authors 



^In the ground state of the even-n clusters with chain geometry, the highest occupied orbital is 
doubly degenerate and we assumed each spatial orbital to be occupied by one electron to make a 
triplet ground state. For other cases, all the occupied orbitals are completely filled. 
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in ref. p8| , p!9| observe transitions in heavier even-n clusters which fall within the same sys- 
tematics as the weak odd-ra transitions. We believe that the even-ra transitions may be 
associated with the single-particle transition across the energy gap. There are two inequiva- 
lent electrons at the Fermi surface due to the partial orbital occupancy, and one combination 
becomes the high-frequency collective state and the other combination remains near the gap 
energy with only a small oscillator strength. In Fig. 4 we show the systematics of the low 
transition, compared to the calculated tt — vr* gap energy. 

We have also calculated the in-plane response of ring configurations in the range = 
7 — 15. For simplicity we assumed uniform circular rings, although the structures may have 
some distortion [^]. The results for the collective tt transition are shown in Table II. The 
ratio of collective frequencies for rings to chains is plotted in Fig. 5. We see that the rings 
are predicted to have a collective frequency about twice that of the chains with the same 
number of atoms. 

III. INTERPRETATION 

Our goal in this section is to characterize the n- dependence of the excitations, and further 
the dependence on shape, whether chain or ring. It is well known that the Fermi gap (or 
HOMO-LUMO gap) Ae in a linear chain depends on n as 

1 

Ae oc -. (6) 

n 

This function fits well the systematics of our calculated single-particle transition energy, but 
it does not describe the collective excitation. This may be seen from the log-log plot of the 
energies in Fig. 6. The figure shows that electron-electron interaction plays an important 
role not only in the absolute frequency of the transition but in its functional form as well. 
The inadequacy of eq. (6) to describe the systematics of the strong excitation was also 
mentioned in ref. . 



In purely classical physics, the analogous problem is the plasma resonance in a needle- 
shaped conductor. The closest problem that can be treated analytically is plasma resonance 
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in a ellipsoidal conductor. Ref. |^T|] derives the formula 

1 - / ^ 1 . 1 + e 



,2 



= — + — log- ]ujI (7) 



e2 V " ' 2e 1 - e 
where e is related to the ratio of short to long axes, R±/R\\, by 



This formula has been applied to the collective vr excitations in football-shaped fuUerenes 
p^ . For our problem, we obtain the fit shown in Fig. 6, treating e and c^o as adjustable 
parameters. It describes the systematics very well. Of course, the carbon chain is not 
ellipsoidal, and we should seek rather that analytic behavior of the resonance in a conducting 
cylinder. 

The behavior of electromagnetic resonances on infinitely long wires is known from clas- 
sical electromagnetic theory. The dispersion formula for the one-dimensional plasmon on a 



long wire reduces to the following expression in the long-wave length, thin wire limit P3| , |24 



2 Anriee^ 2, 1 /n^ 

cu = q log — (9) 

m qa 

where q is the reduced wave number of the plasmon, Ue is the density of electrons per unit 
length, and a is the radius of the wire. For a finite wire, the lowest mode would have a q 
varying inversely with the length of the wire L. Thus the lowest mode would behave as 



uj oc - — . (10) 

This in fact is the asymptotic behavior of eq. (|^) in the limit of large R\\ = L. Taking 
L cc n, the frequency dependence in chains would be 



UJ oc (11) 

n 

This behavior can be extracted from a more quantum approach using the polarizability 
estimate of the collective frequency ||25|| , 

= . (12) 

ma 
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Here N is the number of active electrons and a is the polarizability. This formula is derived 
from the ratio of sum rules, and N may be identified with the oscillator strength / associated 
with the transition. We established in Sect. 2 that the oscillator strength in the tt manifold 
of states is given roughly by the number of tt electrons, and scales accordingly with n. 

We next estimate the polarizabihty. The asymptotic behavior for large n can be deter- 
mined under the assumption that the chain behaves as a perfect conductor. The electrons 
in a perfect conductor will respond to an external field to restore the internal electric field 
to zero. This gives an implicit equation for the electron density in terms of the external 
electric field S, 

z£ = ej dz'6n{z')V{z - z') (13) 

Here 5n{z) is the induced linear electron density, obtained by integrating the induced or- 
dinary density over transverse coordinates. If this is solved for 5n{z), the polarizability is 
then obtained as the ratio of the induced dipole moment to 



a 



= ^ 1 dz'z'Sn{z') (14) 



to determine asymptotic behavior, we note that V{z) is strongly peaked at zero, so we may 
approximate it as a 5 function times / Vdz: 

J dz'Sn{z')V{z - z') ^ 6n{z) J dz'V{z') (15) 

The V{z) behaves as l/\z — z'\ at large separation, so integral on the right hand side depends 
logarithmically on the integration limits. Thus we may estimate the integral as 



/ dz'V(z') / dz'V(z') f« 21og(-). (16) 
J J-L/2 a 



Here a is a length having the order of magnitude of the transverse dimension of the wire. 
We thus obtain 

^ ' 21og(L/a) ^ ' 

10 



From eq. (4), this implies that the polarizabihty is 

" = 241og(L/a) 

We now insert this in eq. (12) and use n ~ L to obtain eq. (pi]). The fit with this function 
is shown in Fig. 7. It obviously describes the asymptotic behavior much better than eq. 
(6), which is also shown in the figure. 

We finally discuss the relative frequencies of the modes in chains and rings using eq. 
(12). We first examine the oscillators strengths, needed for the numerator of eq. (12). 
Naively we would expect similar values, since the number of tt electrons in a ring from 
orbital counting is given by = 2n. However, comparing Tables I and II, we may see that 
the in-plane ring values are only 2/3 the chain values along the axesQ. Much more important 
for the frequency change is the differing polarizabilities of chains and rings. We have made 
a similar asymptotic estimate of the ring polarizabihty, and we find that it also increases 
logarithmically with the circumference of the ring. The ratio of polarizabilties of a thin wire 
ring compared to a straight wire of the same length is about a factor of 6. The two factors 
of the oscillator strength ratio and the polarizabihty ratio combine in eq. (12) to produce 
a frequency shift by a factor of two. This is indeed what the TDLDA gives asymptotically 
(see Fig. 5), confirming the polarizabihty and oscillator strength argument we made here. 

IV. CONCLUSION 

We have demonstrated that the collective vr transition in carbon chains and rings behaves 
in a very systematic way, calculable to ~ 10% accuracy by TDLDA, and understandable in 
macroscopic terms. This should give one confidence in using these transitions to infer the 



^The reason for tlie lower value is that tlie electric field in the plane of the ring is partly transverse 
and partly longitudinal with respect to the C-C bond axes. This implies dipole excitations lie partly 
outside the tt manifold of states. Thus part of the sum rule is shifted to higher energy excitations. 
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shape of the clusters. We hope in addition that the smooth dependence of the transition 
would allow its perturbation going between similar structures to be calculable to accuracy 
of interest for spectroscopic identification purposes. 

We thank R.A. Broglia for calling our attention to eq. (|^. This work is supported by 
the Department of Energy under Grant No. DE-FG06-90ER40561, and by a Grant in Aid 
for Scientific Research (No. 08740197) of the Ministry of Education, Science and Culture 
(Japan). Numerical calculations were performed on the FACOM VPP-500 supercomputer 
in RIKEN and the Institute for Solid State Physics, University of Tokyo. 



APPENDIX: NUMERICAL ASPECTS 

The method we use, explicit time integration of eq. (1) in the coordinate space represen- 
tation, can be shown to be the most efficient of the methods in use for large systems when 
there are no symmetries to reduce the dimensionality. The coordinate space representation 



is also the most efficient for the static problem under similar conditions . 

Let us first consider the matrix RPA in a particle-hole configuration space representa- 
tion. This method has been applied to sodium clusters in ref. and to Cgo in ref. [^. The 



number of numerical operations to extract the eigenmodes and eigenfrequencies scales with 
the dimensionality of the matrix M as M^. The dimensionality of the matrix is given by the 
number of particle-hole configurations. A complete calculation, guaranteed to respect the 
sum rules and conservation laws, requires that all occupied and unoccupied orbitals be in- 
cluded in the space. The single-particle space will have a dimensionality that is proportional 
to the size of the system, and the number of unoccupied orbitals will thus be proportional to 
the number of particles N . Thus the dimensionality with this method scales with as N^. 
Note that this could be much reduced by truncating the particle-hole basis. For example, 
in our study of collective vr transition in carbons chains, the orbitals away from the Fermi 
energy are not so important. However, we do not know a systematic scheme for truncation 
that would preserve the oscillator strength and lead to a more favorable A^ dependence for 
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the algorithm. 

Another widely used technique is the linear response method. This has been applied 
to the molecule N2 in ref. and to alkali metal clusters in ref. [^]. Here the object one 
calculates is the density-density response function. In a coordinate space representation, its 
dimensionality is the number of points in the coordinate space mesh. Thus the dimension- 
ality of the matrix scales as M oc A^. The matrix operation required in this case is inversion 
rather than diagonalization, but it also goes as the cube of the dimension. So this method 
scales as 0{N^). Note that we could have used another representation of the response, such 
as momentum space, and still obtained 0{N^) scaling. 

The advantage of direct coordinate space methods is that one can take advantage of the 
sparse character of the Hamiltonian in that representation. The Hamiltonian has a dimen- 
sionality that scales M oc N, and the number of operations required to apply the Hamilto- 
nian to a single-particle wave function is also 0{N) because of the sparseness. Performing 
the operation on all particles is then 0{N'^), and this gives the scaling for the method. 
Note that there is a large prefactor, because the equations have to be integrated many time 
steps. However, the time integration depends only on the energy scales in the Hamiltonian 
which does not change with A^. Finally, we note that the advantage of this method would 
apply to any representation of the wave function that allows a sparse Hamiltonian matrix. 

For the numerical aspects of integrating the Kohn-Sham equation, we follow closely 



the method used in ref. ||T2| to solve the time-dependence mean field equations in nuclear 
physics. The algorithm must insure energy conservation and norm conservation to very high 
accuracy to be useful. In addition, for the results to be converged, the time step of the 
integration must be small compared to the inverse frequencies of the density oscillations. 
Fixing the time step by this criterion, energy conservation is achieved by using an implicit 
method to integrate the equations. One can show that the integration over the time step 
will conserve energy if the mean field is computed with a fixed density given by the average 
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of the initial and final densities, for Hamiltonians with ordinary two-body interactions^. In 
our computer program, we obtain sufficient accuracy for our time steps by using a simple 
predictor corrector method to find the average density and integrate over the time step. 

We have now reduced the problem to the time integration of a fixed Hamiltonian. Since 
the Schroedinger equation is linear, we can integrate straightforwardly by using the Taylor 
series expansion of the time evolution operator, 

^max 1 

<P{t + At) = e-'^^V(t) = ^ -{iHAtY(t){t) (Al) 

n=0 ^• 

We found that the fourth-order approximant, kmax = 4, is sufficient with our time steps to 
insure norm conservation to the needed accuracy. 



more general condition, valid for nonlinear density dependencies such as in the correlation- 
exchange potential, is given by taking the potential in the single-particle Hamiltonian as 

y^^-^ ^ v[n+ir)] - v[n^{r)] 
— 71- 

where v is the potential energy functional of density and rij- are the densities at the beginning and 
end of the time step. 
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TABLES 



TABLE I. Excitation energy and strength of carbon chains in TDLDA 



Size 


7~i r T rl 

-Efrcc [cV] 


/free 


TP r T /"i 
-C/TDLDAieV/J 


/tdlda 


3 


4.15 


3.5 


8.1 


3.1 


4 


3.89 


5.2 


7.2 


4.5 


5 


2.81 


6.9 


6.4 


6.3 


6 


2.72 


8.6 


5.9 


8.0 


7 


2.13 


10.2 


5.3 


9.8 


8 


2.10 


11.9 


5.0 


11.4 


9 


1.71 


13.5 


4.6 


13.1 


10 


1.71 


15.2 


4.4 


14.8 


11 


1.43 


16.7 


4.1 


16.4 


12 


1.44 


18.4 


3.9 


18.1 


13 


1.23 


19.8 


3.7 


19.7 


14 


1.26 


21.6 


3.5 


21.3 


15 


1.08 


22.9 


3.3 


22.9 


20 


0.90 


31.7 


2.7 


30.8 
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TABLE II. Excitation energy and strength of carbon rings in TDLDA 



Size 


^free [eV] 


/free 


-E^TDLDA [eV] 


/tdlda 


7 


4.9 


5.8 


7.8 


3.1 


8 


4.8 


8.6 


8.2 


4.6 


9 


4.1 


9.0 


7.6 


5.5 


10 


3.7 


9.8 


7.2 


6.9 


11 


3.5 


10.7 


7.0 


7.9 


12 


3.3 


11.8 


6.9 


9.2 


13 


3.0 


12.5 


6.6 


10.3 


14 


2.8 


13.5 


6.3 


11.3 


15 


2.6 


14.4 


6.1 


12.4 


20 


2.0 


18.7 


5.2 


17.4 



18 



Figure Captions 

Fig. 1. Optical response of the linear chain Cy, calculated in TDLDA. The mesh spacing 
is Ax = 0.3 A and At = 0.001 eV^^. The spatial grid has the shape of a cylinder, with 
30,000 mesh points. The integration time was 10 eV~^. 

Fig. 2. Molecular orbital diagram for the Cy chain in the LDA. On the left are shown the 
lowest IT states, which correspond to the 7 orbitals of the tight-binding approximation. On 
the right are the highest occupied and lowest unoccupied a orbitals. 

Fig. 3. Systematics of the longitudinal mode in Cat linear chains. Circles: experimental, 
from ref . fl^ and ; squares, TDLDA; triangles, quantum chemistry calculations from ref . 
g and Inl. 



Fig. 4. Systematics of vr — vr* gap energies compared with observed transitions from ref. 
11,111 • 



Fig. 5. Ratio of the collective excitation energy in rings with compared to chains, from the 
TDLDA. 

Fig. 6. Systematics of the single-electron transition and the collective excitation in chains, 
compared with functional fits of the form eq. (6) and eq. (7). 

Fig. 7. Comparison of the TDLDA collective excitation in chains with eq. (6) and (11). 
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